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Abstract. The variety of magnetic phases observed in rare-earth heterostructures at low 
temperatures £Q, such as Ho/Y, may be elucidated by an ANNNI-like model Hamiltonian. 
In previous work modelling bulk Ho 0, such a Hamiltonian with a one-dimensional 
parameter space produced a single multiphase point. In contrast, the parameter space 
of the heterostructure model is three-dimensional, and instead of an isolated multiphase 
point, we find two-dimensional multiphase regions. In an example of Villain's "order from 
disorder" [3J E] , an infinitesimal temperature breaks the ground-state degeneracy. In first 
order of a low-temperature expansion, we find that the degeneracy is broken everywhere in 
a multiphase region except on a line. A segment of the line appears to remain multiphase to 
all orders in a low-temperature expansion when the number L of magnetic layers between 
non-magnetic spacers is 4 but not for other values of L. For L = 4, the hierarchy of phases 
more closely resembles that in the ANNNI model than in the bulk six-state clock model on 
which the present model is based. 
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1. Introduction 

Layered planes of rare-earth metals exhibit a wealth of magnetically-ordered phases at low 
temperature. In helimagnetic phases, spins (treated classically) align ferromagnetically 
within each plane, with an axial RKKY interaction responsible for a progression of spin 
angles through successive planes Strong easy-axis anisotropy may frustrate the natural 

RKKY pitch angle, leading to a multitude of possible phases characterized by the number of 
layers separating skips, or "walls," in the pattern of pitch angles. In the axial-next-nearest- 
neighbour Ising (ANNNI) 013101110] and related clock models PH Q21 EEH1 El El EES] , a single 
parameter controls the relative strengths of competing interactions, and at a single value 
of this parameter, infinitely many phases coexist; this is called a multiphase point. Since 

f Current address: Department of Mathematics, University of Florida, Gainesville, FL, USA 
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these phases cover all allowed spacings between walls, such phases are indistinguishable from 
random sequences. Thus the zero-temperature state is disordered. This disorder is broken 
at infinitesimal temperature in an example of "order from disorder" [HI E] . We now ask what 
happens in a model of helimagnetic heterostructures with a three-dimensional parameter 
space: we identify fully two-dimensional multiphase regions and investigate the topology of 
the low-temperature phase diagram. 

With the giant magnetoresistive effect JH] in ferromagnetic/nonmagnetic superlattices 
having spawned important technological applications that reached the market around 1997 
fTJEI], it seems practical, as well as theoretically interesting, to examine the possible phases 
of helimagnetic/nonmagnetic superlattices. Such superlattices have been deposited using 
molecular-beam epitaxy, alternating dysprosium erbium [20], or holmium pQ with non- 
magnetic yttrium spacer layers as well as holmium with lutetium j2l] . Surprisingly, neutron- 
scattering experiments show that the helicity of the spins in the rare-earth layers is preserved 
across the spacers, with the magnetic moments forming long-period "spin-slip" phases pQ. 
RKKY-like polarization [22] of conduction electrons in the non-magnetic layers is again 
implicated [2ID 121]; m any case, we can model the indirect exchange across non- magnetic 
spacers in parallel with that between successive magnetic planes. If the exchange parameters 
can be controlled with pressure, external fields, or spacer-layer thickness, such systems could 
possibly be useful as magnetic sensors or in data-storage applications. Axially modulated, 
high-order, commensurate phases are not limited to rare-earth heterostructures: Szpilka and 
Fisher jS] cite half a dozen other systems in which such phases have been observed, ranging 
from CeSb j2S] to ferroelectric thiourea [2EII2Z!- 

Seno et al. [2 j applied the ANNNI ideas to a case of infinite hexagonal anisotropy, the 
six-state clock model, relevant, for example, to bulk holmium.f A spin a in the j th plane 
points in a direction that is an integral multiple, rij a , of 27r/6. At zero temperature, all the 
spins in a plane point in the same direction (rij), and the model is controlled by a single 
parameter, the ratio x of the strength of the next-nearest-axial-neighbour antiferromagnetic 
(J2) to nearest-axial-neighbour ferromagnetic (Ji) interaction, with the axial terms in the 
Hamiltonian summing — J\ cos(27r(n J+ i )Q , — n_,- ja )/6) and + J2 cos(27r(n :) ' + 2 ) a — 7ij jCe )/6). For 
< x < 1/3, the ground state is a ferromagnet, for 1/3 < x < 1 a helimagnet with 
no walls, and for x > 1 a helimagnet interrupted by walls every second layer. At the 
single point x = 1 in the one- dimensional phase diagram, infinitely many phases coexist 
in the ground state. We represent the helimagnetic phase (1/3 < x < 1) by the axial 
sequence . . . 012345012 . . . , understanding that this includes as well the translations and 
reflections of the sequence. The two coexisting period-2 phases for x > 1 are represented 
by . . . 00330033 . . . and . . . 01 1 34 1 01 1 34 . . .: this last is thought of as a modification of the 
helical phase by the insertion of skips, or walls (denoted "|"), every second layer. The walls 
are analogous to domain walls in the ANNNI model. At the multiphase point, x = 1, in 
addition to . . . 00330033 . . . , a helical phase with walls placed anywhere at least two layers 
apart is a ground state of the system, e.g., ... 01 1 345 1 12 1 450 .... A convenient notation in 

f An extension of this work presented a small-inverse-anisotropy expansion about the clock model and again 
found a hierarchy of phases emanating from the multiphase point at infinite anisotropy |28| . 
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ANNNI-type models labels a periodic phase by the spacings between successive walls: thus, 
this last example is (23), the phase with walls every second layer (2), and the bare helical 
phase without walls (oo). In a low-temperature expansion, Seno et al. followed a hierarchy 
of phases (similar to what we describe below) and showed that each phase between (23) and 
(oo) acquires a region of stability at infinitesimal temperature. 

The forgoing model simplifies the actual magnetic structure of bulk holmium. Neutron 
scattering gives the turn angle per atomic layer as 30° rather than 60°, with moments bunched 
in pairs around the six easy axes (2H1 U\, and while the average turn angle increases in 
films, the effect is thought to be due to interspersal of singlets among the pairs; thus the 
(3) phase in the simplified model might actually represent moments . . . 00122344. . ., where 
pairs of repeated spins lie a few degrees before and after the easy-axis direction (see Fig. 14 
of Reference P). The model, or its present extension to superlattices, was meant not to 
reproduce realistic details of a particular rare-earth helimagnet but rather to reduce a system 
with competing crystal-field and exchange interactions to the simplest form, in which exact 
results are possible, so as to investigate universal properties of the resulting hierarchy of 
commensurate, longitudinally-modulated spin-slip phases. 

2. The model and its ground states 



We consider a superlattice in which blocks of L magnetic layers are separated by non- 
magnetic spacers characterized by effective couplings J[ and J 2 , this simple extension of the 
bulk model of [2] gives the full Hamiltonian 




where i labels layers, a a spin within a (simple-hexagonal) layer, and /3(a) its nearest 
neighbours. The unprimed sums in the second line are taken only over bonds that do not 
straddle a non-magnetic spacer, while the primed sums in the third line are taken only over 
bonds that do. For purposes of the low-temperature expansion, the in-plane ferromagnetic 
coupling constant Jo is taken to be positive and much stronger than any of the axial couplings 

Since we are looking for helical phases, we take all of the remaining four couplings 
also to be positive. (Certain negative couplings are in fact related to the positive sector 
by symmetries of TC.) The model reduces to that of [2j when J[ = J\ and J' 2 = J2 or, 
equivalently, when L — 1. The three-dimensional coupling space is given by x = J2/J1, 
y = J[/Ji, and z = J' 2 / J\\ it is convenient to set J\ = 1. 

We generalize the previous notation to accommodate states of a superstructure in which 
blocks of L magnetic layers are separated by non-magnetic spacers, denoted by || , with the 
arrangement repeated periodically. (The symbol || may denote any number of atomic layers 
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of the non-magnetic metal.) Since the direct interactions in extend a maximum of two 
layers in the axial direction, walls are classified in three categories. A wall at least two layers 
from a spacer has the same energy cost as in the bulk model and is termed a type-1 wall, 
for example (L = 5) 

... || 0123 1 50 1| 12... . (2) 

Insertion of a wall one layer from a non-magnetic spacer has a different energy cost, since a 
J' 2 bond is broken. This is termed a type-2 wall: 

... || 01234 1 1| 12... . (3) 

A type-3 wall coincides with a non-magnetic spacer: 

... || 012345 1| 1 12... . (4) 

Helical configurations, including (oo) itself, that differ from (oo) only by the insertion of 
walls are called wall states. These states preserve the sense of helicity (positive or negative). 
We consider L > 3, as L = 1 is the same as bulk, while L = 2 omits the J2 (x) parameter and 
so has only a two-dimensional parameter space. It is also less likely to be of experimental 
interest. 

A straightforward calculation yields the total energy of a wall state as a function of the 
densities Wi of walls of the the three types: 

~ ((1 + x){L - 2) + 1 + y + 2z) + (1 - x)Wx + (l - ^) W 2 + (y - z)W 3 . ^ 

As in the original model, successive walls are energetically forbidden. We seek regions of the 
three-dimensional parameter space in which the insertion of a wall of some type costs no 
energy: this occurs when the coefficient of one or more of the densities Wi vanishes. Thus 
the planes x = 1, (x + z)/2 = 1, and y = z all potentially constitute multiphase regions; 
however, it is also necessary to consider competing non-wall states, which may have lower 
energies. For present purposes, we shall concentrate on the y = z plane, for which type-3 
walls cost no energy. Since a negative energy for type-2 walls would shut type-3 walls out, we 
examine the part of the y = z plane to the left of the x + z = 2 line. By considering points to 
the left of the line x — 1, we exclude type-1 walls as well. For L — 4, direct calculation gives 
the phase diagram of Figure H An exhaustive computer search (of phases of length 3L = 12 
with twisted periodic boundary conditions) verified that the wall-state energy © is lower 
than that of any competing phase inside a triangle in the y = z plane, which constitutes a 
multiphase region. Comparable wall-state regions were calculated numerically for L in the 
range 3-11. 

Since we are concentrating on a region in which type-1 and type-2 walls are excluded, 
while type-3 walls cost no energy, we adapt the notation of [2] to count magnetic blocks, 
rather than magnetic layers, between walls. Thus, for example, with L — 4, (1) has a wall 
coinciding with each spacer, while the (2) phase has a wall at every other spacer. Since no 
restriction prevents adjacent walls of this type, the count of possible phases is simply 2 to 
the power of the number of magnetic blocks; this represents a simplification relative to the 
ANNNI and other related models [SOJ - 
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Figure 1. Ground-state phase diagram for (part of) the y — z plane, L = 4. The horizontal 
axis gives the normalized bulk second-neighbour coupling, the vertical the couplings across 
non-magnetic spacers. Outside the triangle delimited by dot-dashed lines, the ground states 
are as indicated. Inside the triangle, wall states are the ground states. (An exhaustive search 
found no lower-energy states of length up to 3L.) Similar ground-state phase diagrams were 
calculated for other values of L. The first-order low-temperature expansion gives (oo) to the 
left of the dotted line within the triangle and the (1) phase to the right; on the line itself, 
these phases and their progeny coexist, requiring a higher-order low-temperature expansion 
to distinguish. On roughly the upper half of the left leg of the dotted line, from z = 11/9 
to z — 13/9, we believe infinitely many phases of the form (l fe 2) coexist to all orders. 

3. Low-temperature expansion: first order 

The novel feature presented by the current problem is the multiphase triangle (for L = 4 
or a similar polygon for other L) throughout which infinitely many phases coexist at zero 
temperature. An interesting theoretical question is how thermal disorder can distinguish the 
free energies of all these phases in the given region. 

Although the Hamiltonian contains only first- and second-neighbour axial terms, a 
non-zero temperature introduces effective long-range interactions through an axial chain of 
thermally-excited spins, each pointing in a direction at variance with its in-plane neighbours 
El Ej By analogy to the ANNNI model, we call such excitations "spin flips." Since 
the number of ways an excitation of a particular energy may occur depends on the state, 
flipped spins provide an entropic mechanism for distinguishing the free energies of wall states 
at infinitesimal temperature. If the i th excitation, which may involve several spins, has an 
energy AEi relative to the ground-state energy per spin E and can be placed on the lattice 
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of N spins gi different ways, the free energy per spin is given by the linked-cluster theorem 
jSU: 

f = E -k B Tj2^- pAEl , (6) 

i 

where ji = Ubin^o g%/N is the intensive part of gi/N. (The limit discards those terms 
in gi that go as higher powers of N; such terms come from independent clusters of spin 
excitations.) 

We apply the method first to an isolated spin flip, which may occur in a layer adjacent 
to or one layer separated from a spacer, or it may (L > 4) occur in bulk. An isolated spin 
flip in bulk gives the same contribution to / regardless of phase, so we calculate the energies 
and counts 7$ just for the first two cases, leading to the weights in Table H The case L = 3 
requires special treatment because the cost of an excitation in the layer in the middle of a 
block depends on the presence or absence of walls on both sides. 

Table 1. Contributions to JfjJ) are formed by a count (per spin) of the number of ways of 
forming the excitation times a Boltzmann factor. The left column gives an example of the 
excitation under consideration, where the caret (A) marks the plane in which a single spin is 
rotated ("flipped") plus 60° or minus 60° from the angle of its neighbours in the plane. The 
second column gives the Boltzmann factor, and the remaining columns give the intensive 
counts 7j weighting the Boltzmann factor for the cases (1), (2), and (oo). L is the number 
of magnetic layers in a block. The last three rows apply only to L = 3. Here, (3 is the 
inverse temperature, q = exp(— @Jq/2), t[— 6] the number of in-plane nearest neighbours, 
and r = exp(— (3J\j e 2). 
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We consider L > 4 first. If there are no type-3 walls, the only single-spin excitations 
(other than bulk) will be of one of the types in the first two rows of Table ^ This describes 
the (oo) phase. If a phase has the maximum density of type-3 walls, the excitations will be 
of the types in the second two rows. This is the (1) phase. To this first order in the low- 
temperature expansion, any other wall phase (e.g., (2)) will have a free energy intermediate 
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between these two cases. Thus we look first for the coexistence of (1) and (oo). Subtracting 
rows 1 and 2 from the sum of rows 3 and 4 gives the free-energy difference 

= f(l) - /(oo) 

9 

ksTq 1 (r l ~ x+z + r 1+2x+z _|_ r 2-x-2y+z _)_ r -i+2x+y+z ^ 

L 

_ r l-x+2z _ r l+2x~z _ r .2-x-y+2z _ r ~l+2w+2y-z\ 

where (3 = 1/(&bT) is the inverse temperature, t the number of in-plane nearest neighbours, 
q = exp(— (3J /2), and r = exp(—j3Ji/2). Setting Af = and y = z yields the expression 

r z +r _i_ r l-2z 

r = . (8) 

In the zero-temperature limit, r — ► 0, so the power of r with the smallest exponent dominates. 
This allow us to solve for the coexistence line, 
2 

for < z < ~ 
3 1 

l--z for | < z < 1 (L > 4), (9) 



x 



3 

1 



2 

for 1< z < 2 



3 

drawn as a dotted line in Figure H In the multiphase region to the left of this line, the 
(oo) phase has the lowest free energy, breaking the infinite degeneracy of zero temperature. t 
To the right of the line, the (1) phase dominates. On the line itself, all wall phases remain 
degenerate; to break the degeneracy it will be necessary to consider more flipped spins. 

First, however, the model with L = 3 introduces a new element to the low-temperature 
expansion. In the last three rows of Table ^ the count of the (2) phase does not merely 
interpolate between the counts of (1) and (oo): that is, a single-spin excitation in the middle 
plane of a magnetic block distinguishes not only (1) from (oo) but also each from (2). Thus, 
the first-order expansion must potentially consider three coexistence lines. In the event, the 
three collapse to one. For z > 0, all wall phases coexist on the line 

z=\{l-x) (L = 3). (10) 

For z > (3/2)(l — x), the (1) phase has the lowest free energy, while for smaller z, the (oo) 
phase has the lowest free energy. 



4. Expansion to higher orders 

The hierarchy of potential phases in the low-temperature expansion has been described well 
elsewhere El El IH2] and so will only be summarized. At any order of the expansion, a 
coexistence region has been established between two "parent" phases and infinitely many 
other wall states. (In Figure ^ for L — 4, this region is the zig-zag line, on which, to first 
order, parents (1) and (oo) coexist with all other wall states.) Spin excitations to this order of 

f The twelvefold degeneracy of (oo) neither scales with iV nor affects the spin-spin correlation function. 
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the expansion cannot distinguish the parents from the other wall states, but by adding some 
number of additional spin excitations, linked to those of the given order, we can distinguish 
the two parents from a "child" phase made by concatenating one period of each parent. As 
examples, the child of (1) and (oo) is (2), while that of (1) and (2) is (12). A connected 
chain of spin excitations can "see" the presence or absence of walls over its length; viewed 
another way, this leads to an effective long-range interaction between walls. 

While in principle one could continue the enumeration of connected excitations of two, 
three, and more spins along the lines of Table [""J a transfer- matrix technique ["QJ E] is well 
suited to computer symbolic algebra. We defer implementation details to the appendix. The 
matrices are more involved than those in so the results are for specific cases, from which 
we conjecture generalizations. 

In first order, we have already seen the two-dimensional multiphase region shrink to one 
dimension (Figure EJ. We wish to find out whether the line shrinks further to a point or set 
of points, or whether the line, or a portion of the line, behaves like a multiphase point, with 
the additional degree of freedom essentially irrelevant. It is also of interest whether all wall 
states descending from (1) and (oo) attain stability or only a subset. 

We carried out the low-temperature expansion for magnetic blocks of length L between 
3 and 17; except for the interesting case of L = 4, the hierarchy terminates after just a few 
phases. Aside from (oo), the only stable phases found for L = 4 were of the form (l k 2), 
< k < 27 (the highest calculated) and k = oo (i.e., (1)). This resembles the ANNNI 
model [""J IB] more than some clock models in that there do not exist two phasesf all of whose 
progeny attain stability. Villain and Gordon [SSj (see also [""]) distinguish a Devil's staircase 
[3*4] from a "harmless" one. In both, a multiphase point gives rise to a large number of 
phases that approaches infinity at T — > 0. However, in the latter case, at any finite T > 0, it 
is argued that only finitely many phases are stable. Since our model fails to find an infinite 
hierarchy of "mixed phases" [12], we conjecture that our staircase may similarly be harmless. 

The way the (l)-(oo) coexistence line breaks up for L = 4 is also of interest. It intersects 
the multiphase triangle (Figure H} for 1/3 < z < 13/9; outside this region, it ceases to 
describe coexistence of ground states. The symbolic transfer-matrix calculation finds that 
(2) is stable on the line only for 3/4 < z < 13/9. Below 3/4, there is a first-order phase 
transition between (1) and (oo). The phase (12) is stable at z = 3/4 and then again for 
11/9 < z < 13/9. All subsequent phases (l fc 2) for which we were able to extract symbolic 
results (k < 5) are stable for 11/9 < z < 13/9 (that is, 3/4 and 11/9 drop out). Numerically, 
k =6-27 is stable for 11/9 < z < 13/9, the "<" indicating the inability of the numerical 
code to distinguish between the proper and improper inequality. We conjecture that (l fc 2) is 
stable for 11/9 < z < 13/9 for all k > 2 and that "mixed phases" never come in, something 
we were able to confirm up to the mixed phase (l 18 2 l 17 2). 

For L = 3, the coexistence line (|1U|) intersects with the region in which wall states have 
the lowest energy for < z < 1. The states (1), (oo), (2), (12), and (3) are stable on this 
line segment, but no other phases. 

f (23) and (oo) in the bulk model [5] 
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For L = 5, the coexistence line is again ®, which passes through the wall-state region 
for < z < 13/9. The same phases are stable as for L = 3: (2) for < z < 3/4, (12) for 
< z < 3/4, and (3) for < z < 3/4. For L = 6, the (l)-(oo) coexistence line gives the 
lowest energy for < z < 13/9; however, the only other stable phase is (2), and only at the 
single point z = 3/4. 

The following pattern appears to hold for L > 6: the coexistence line (JSJ) intersects with 
the wall-state region for < z < 13/9. Even values for L (we computed 8, 10, 12, 14, and 
16) give a first-order transition between (1) and (oo) all along the coexistence line. No other 
phases are stable. For odd L (7, 9, 11, 13, 15, and 17), the phases (2), (12), and (3) are also 
stable for < z < 3/4. 

5. Implications 

The low-temperature expansion applies at infinitesimal temperature, but the bulk model 
has also been investigated with a mean-field theory, which should be valid only at high 
temperature |2j. The low-order phases predicted by the low-temperature expansion were 
seen to spread out from the multiphase point as temperature increased (Figure |2J); the only 
notable discrepancy between the two extreme theories was the presence of phases (2 fc 3) in 
the mean-field calculation, and this was explained in terms of a competing phase. Near 
the zero-temperature multiphase point, which in the L = 4 model would be replaced by 
the multiphase zig-zag line of Figure Q the spin-spin correlation length is expected to be 
small, as a large number of commensurate phases with different periods coexist. At any 
temperature greater than zero and less than the Curie temperature, only one phase is stable; 
however, in an experimental system, interfacial roughness and interdiffusion might lead to 
coexisting commensurate phases from nearby points in the temperature-phase space. As the 
temperature increases, the volumes of stability do as well, so that no phases lie nearby, thus 
stabilizing a single phase. 

Interestingly, the coherence lengths £ of the basal-plane holmium moments in Ho/Er 
superlattices have been found to increase with temperature T between 8K to 100K [35J. 
Since Er acquires a moment below 100K, the experimental system is considerably more 
complex than our simple model; moreover, too few temperatures were measured to permit 
a comparison to the plateaux one would expect in £(T) from Figure El A similar effect is 
observed in Er/Lu [3BJI2H]- 

The question of commensurate versus incommensurate magnetic modulation also awaits 
experimental resolution. In the low-temperature expansion, incommensurate phases are 
only approached, as the limit of a hierarchy of commensurate phases, while the bulk mean- 
field calculation (Figure |2j) suggests that these limiting phases will occupy a volume of 
measure zero in the phase diagram. While several rare-earth systems unambiguously show 
commensurate phases [2HJ 12D1 EH] , other superlattices appear to show a continuous increase 
with temperature in the average turn angle per atomic layer, suggesting that incommensurate 
phases are generic |37|. We cannot rule out an averaging effect being responsible, 
but this would appear inconsistent with the absence of plateaux and the expectation of 
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Figure 2. Diagram of a hierarchy of phases emanating from a multiphase point at zero 
temperature up to a Curie temperature. The horizontal axis represents a ratio x of coupling 
strengths, which at 1 leads to zero-temperature disorder. Raising the temperature from the 
vicinity of x = 1 gives a succession of stable phases. Adapted from Figure 1 of where it 
shows a numerical mean-field calculation on the bulk six-state clock model. In the present 
context, it can be thought of as schematic for the T > behaviour of the system of Figure ^ 
at some point along the zig-zag line, where x represents a transverse dimension. 

vanishing measure for high-order phases. It will be particularly interesting to investigate 
whether a statistical-mechanical model not much more complex than that considered here 
can incorporate more of the qualitative behaviour seen in rare-earth superlattices. 

We have shown that a superlattice of helimagnetic and non-magnetic layers exhibits 
behaviour different from that of the bulk six-state clock model |2J. There are multiphase 
regions, rather than a single multiphase point. When precisely four magnetic layers lie 
between non-magnetic spacers, a line segment in the multiphase triangle appears to support 
a set of phases more like that in the ANNNI model E] than like the bulk six-state clock 
model. For other values of L, the low-temperature expansion finds only a few stable phases. 
This raises the interesting experimental question of whether rich magnetic phase diagrams in 
artificial superlattices could appear for certain magic spacings while being absent for others. 
If the phase diagram were to depend as sensitively on L as in our model, it might be difficult 
to grow films sufficiently uniform to test the hypothesis; however, if the extent of the magic 
coupling were broader (say, L=4-6), the effect could be observable. Further, a multiphase 
region of coupling space might be more amenable to experiment than a multiphase point that 
requires exact tuning; such a region, however, would need to have the full dimensionality of 
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the coupling space, something we have not yet constructed. 
Note 

In the published version of this paper, the reference at the end of the first sentence of 
Section El reads [12]. It should be 0. 
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Appendix 

In order to calculate the free-energy difference of a child from its parents, we adapt the 
transfer-matrix technique fTJ QUI |2] to the region with only type-3 walls. We begin in a 
region over which, to the order already calculated in the low-temperature expansion, parent 
phases (a) = (a x a 2 . . .) of period p a = £V and free energy per spin f a and (b) = (bi b 2 ■ ■ ■) 
of period ph = J2i h and free energy per spin /& coexist and have lower free energies than 
their parent phases, f We then seek the double free-energy difference 

a(ab) = f(ab) "7 f(a) -7 f (b) (A.l) 

Pa+Pb Pa + Pb 

to leading order. If a/gb) < 0, the child phase (ab) acquires a region of stability. Isolated 
spin rotations (as in Table cannot determine the sign of (|A.1|) . since the three phases, 
(a), (b), and (ab), have the same free energies to first order. We must consider connected 
spin excitations: in general, the Boltzmann weight of two (or more) spin rotations that share 
an axial bond will differ from the weight of the same rotations situated in their respective 
planes such that they do not share a bond. Since the Jo (in-plane) bond is assumed the 
most expensive to break, the shortest excitation that distinguishes (ab) from its parents 
provides the leading term in the low-temperature expansion. This requires that the connected 
excitation should span (jp a + pt, — 1) blocks of length L, in the sense that bonds on each end 
extend through the terminating spacer layers and so sense whether these spacers coincide 
with walls. The transfer-matrix technique keeps track of all the combinations of connected 
and disconnected excitations of this length. 

As in |2j, two cases arise. When the product (jp a + pb — 1) ■ L is odd, an excitation of 
connected spins every second layer distinguishes the child from the two parents, and 2x2 
matrices suffice. When the product is even, we shall need 4x4 matrices. 

The principles are best illustrated by an example. Consider distinguishing (2) from its 
parents (oo) and (1) when L = 5. In the following diagram showing just over one period 

f The period of (oo), for the purpose of (|A.1J) . is 1. 
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of the (2) phase, S represent a magnetic layer, while S represents a magnetic layer with a 
flipped spin: 

S ||| SSSSS || SSSSS . (A.2) 

In the (2) phase, the two extremal spacers ( || ) coincide with walls. In the (l)-phase parent, 
all three spacers coincide with walls, while in the (oo)-phase parent, there are no walls. The 
pictured connected spin excitations, spanning pm + p^ — 1 = 1 block, is the shortest that 
is possible for (2) but impossible for either parent. 

The energy difference (jA.ljl subtracts from the free energy of diagram ()A.2j) the parent- 
diagram free energies. We accomplish this with a product of vectors (lowercase Greek letters) 
and matrices. For this example, we get 

a {2) oc (/3 f - a ] ) A (a - (3) , (A.3) 

where a represents a diagram SSS \\ S, [3 the diagram SSS \\ \ S, and A the diagram SSSSS. 
The duality operator, defined for vectors by v* = (Qv) T , with Q having —1 all along the 
antidiagonal, describes the reversed diagram, e.g., a 1 = S\\SSS, with the clock directions 
also reversed. Since the spin in an excitation can be rotated 60° counterclockwise (+) or 
clockwise (— ), both conditions must be accounted for. The four entries of a matrix stand 
for the four ways the two connected spins in a matrix diagram can be flipped: 

+ : tx) ■ < A ' 4 > 

The entries of a row vector are (H — ), those of a column vector (~J , so that each contraction 
in a matrix product sums over the possibilities for a single spin. Each 2x2 matrix entry 
gives Boltzmann weights for connected and disconnected combinations of the two constituent 
spins, as illustrated in Figure lATk . 

In addition, each matrix entry is a difference between the connected Boltzmann factor 
and the disconnected factor, as specified by the linked-cluster theorem, ©. Vectors terminate 
the product fFigure [Alb ). The following 2x2 matrices are required; common factors of 
q f are omitted, since only the signs of the matrix products in the zero-temperature limit 
matter. 

_ 2x , i_ £n x { (A.6) 
f 2 2(1 — r ) T22N — r ) I 

SS\\SSS C= r z ( 1 1 _ 22 ' , z } ) (A.7) 

\ 1 — r Az 1 — r z J 

r |(l-)(l_ r -2^) r f(l-^)(l _ r x) J ( A - 8 ) 




Figure Al. A matrix element represents two flipped spins, a vector element one. Boldface 
bonds are counted at full strength in the Boltzmann weights, while each of the other bonds 
is counted in two different diagrams and so comes in at half strength, (a) A 2 x 2 matrix 
represents flipped spins (S) in the second and fourth planes, (b) A (column) 2-vector 
contracts with a 2 x 2 matrix to its left, (c) A 4 x 4 matrix represents a flipped spin in one 
(and only one) of the first two layers and in one (and only one) of the second two. (d) A 
(column) 4- vector contracts with a 4 x 4 matrix to its left. (Adapted from Reference 



/ r §(l-z)(-i _ 2z\ §(l-«V r 3«_ r 2z\\ 

SS\\\SSS E=r*( r J_ r2z } /_ r2z j J (A.9) 



i / r 2z 



r 



i / r 



SSS || S a = r* (A.10) 



SSS\\\S (3 = r-2 (A.ll) 



The following environments occur only when L = 3: 



SHI Skills 



r 2 a (1 



-2x^ 



f 2 



(A.12) 



s\\ \sss\\s 



swsssws 



n 



i 



f 2 ^1 — r x 

3z 

r 



2x\ 



t 2 (1 — r 



T 2 

,,2 



1 — 7" _2a; ) + f (1 — t" 3 ^ 



(A.13) 
(A.14) 



When (p a + pfe — 1) ■ L is even, there is no unique shortest leading-order diagram on 
the model of (jA.2|) . Rather, a family of such diagrams with flipped spins every second 
layer except for one pair of axially adjacent flipped spins all span the requisite distance. To 
account for a single adjacent pair anywhere along the length of an excitation, Seno et al. |2] 
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(A.15) 



introduced 4x4 transfer matrices of the form 

/+0 + +0 - +00+ +00- \ 
-0 + -0 - -00+ -00- 
+ +0 + -0 + 0+ + 0- 
\0 -+0 --0 - 0+ - 0-/ 
each entry of which considers four adjacent planes in which a spin has rotated in the positive 
(+) or negative (— ) clock direction, or not rotated at all (0). See Figure lA"Tb . The four 
entries of the upper-right quadrant contain no connected spin excitations and so vanish. 
End-cap vectors (Figure lAlb O acount for the final pair of planes, one of which will contain 
a spin flip. The following matrices and end-cap vectors result (again, the common factor of 
q l is omitted): 



SSSS A- 



1 r(l — r x ) 


— - 1 3r 

rf> 2 1 ty 














r 2 (1 — ip—2x\ 
r2+ z (l — 


r(l - 
















-1) 


r(l 








^ iy 2aj/y,2 7^) 




— r) 


ri(l 




r(l — r x ) 


J 



(A.16) 



SSS\\S B- 



SS\\SS C-- 



SSS\\\S D- 



SS\\\SS E-- 



/ t 2 2(1 — r^) 
7>3~§ (1 — r~ 2x ) 

^ — Tf 1\ 

f 2 [T — T I 
f 2 













\ 



\ y 2 ' >■>■■ 



r2 



/ ^,22-1^ _ r ^ 



T 2 



+2 



(1 



r 

1) 

r z ) 



i / 2? 

7" 2 I f 

7"^ {^f 



1 . 

3z\ 



r z {l 



r 



-2z\ 



5z 7z 

t 2 — r 2 

3 z 



r 



1) 



,2+22/^.22 



\ r* _1 (l -r 2 ) 
/ r^~ z (\ — r x ) 

y 7- —a; (7"l — f 5 j 



„22\ 



/ r* _1 (l 

y r 22-l^ _ r 2 2 -j 



r 22 _ J_) 



1 

7" 2 



! (1 



,3;r 



7" 2 2 \ T — r 
7" 2 (1 — T x 

2r+^/ — i 
if^^^ 2 if 2 



f 2 






^>2 2 ^ ^ 

22n „,-2z 



r 2 2(1 — r z ) / 

\ 




7- 2 + 2z (1 — T" 2 ) 

^ _ r -22^ ^22-1^ _ r 2^ y 













\ 



T 2 



„22> 



r 2 (1 — 7" 

2 + "2" (t" 2 — 1) f2~ Z (\ 



1) 



r2 (r 2 — 1) 

r 2-22^ 
2-2 (1 
2z 



r 



- r 
r 



2^ 



/•- 1 "~(1 



r 






r^(r~^ 2 — 1) r^^r 2 — 1) 



r 2z )J 
\ 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



SS\\SS a 



,1-2 



,1 + 2 



\ r z 2 / 



(A.21) 
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SS\\\SS b = 



/ f z 2 \ 
J.1+Z 



(A.22) 



A straightforward computer algorithm generates the relevant sequence of matrices; as 
one example, for L — 4, 

a {23) = (- a t + tf)ACAEACA(a - b) . (A.23) 

The programme expands and symbolically determines the leading behaviour of a in the zero- 
temperature limit; if a is negative, the child attains a region of stability with respect to its 
parents. For sufficiently long chains of matrices, it was impractical to expand the matrix 
products, and a numerical approach was substituted. 
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